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ABSTRACT 

We have carried out optical (R band) intraday variability (IDV) monitoring of a sam- 
ple of ten bright low energy peaked blazars (LBLs). Forty photometric observations, 
of an average of ~ 4 hours each, were made between 2008 September and 2009 June 
using two telescopes in India. Measurements with good signal to noise ratios were typi- 
cally obtained within 1-3 minutes, allowing the detection of weak, fast variations using 
N-star differential photometry. We employed both structure function and discrete cor- 
relation function analysis methods to estimate any dominant timescales of variability 
and found that in most of the cases any such timescales were longer than the duration 
of the observation. The calculated duty cycle of IDV in LBLs during our observing 
run is ~52%, which is low compared to many earlier studies; however, the relatively 
short periods for which each source was observed can probably explain this difference. 
We briefly discuss possible emission mechanisms for the observed variability. 
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1 INTRODUCTION 

Blazars are a subclass of radio-loud AGNs characterized 
by strong and rapid flux variability across the entire EM 
spectrum and strong polarization from radio to optical 
wavelengths. Microvariability, or intraday variability (IDV) 
is commonly observed across much of the electromagnetic 
(EM) spectrum of AGNs but is particularly common in 
blazars. A change of flux of ~1— 15% within a few min- 
utes to hours reflect extreme physical conditions embed- 
ded in small sub-parsec scales. According to the usually ac- 
cepted orientation based unified model of radio- loud AGNs, 
blazar jets usually make an angle ^ 10° to our line-of-sight 
l|Urrv fc Padovarnl Il995l ). The Doppler boosting of the jet 
emission means that most of what we see from blazars arises 
in those jets. These jet dominated AGNs provide a natural 
laboratory to study the mechanisms of energy extraction 
from the vicinity of central supermassive black holes, the 
physical properties of jets and perhaps also accretion disks. 

The radiation of blazars across the whole EM spec- 
trum is predominantly non-thermal. At lower frequencies 
(through the UV or X-ray bands) the mechanism is almost 
certainly synchrotron emission while at higher frequencies it 
is ve ry probably dominated by Inverse-Compton (IC) emis- 
sion (sikora fc Madeiskill200ll : lKrawczvnskil2004 ). The spec- 
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tral energy distri butions (SEDs) of blazars have a double - 
peaked structure l|Fossati et al.lfl99 8; Ghisel hnT et aUll997h . 
Based on the location of the peak of their SEDs, blazars 
are often sub-classified into the low energy peaked blazars 
(LBLs) and high energy peaked blazars (HBLs); the first 
component peaks in the near-infrared (NIR) /optical in case 
of LBLs and in the UV/X-rays in HBLs, while the second 
component usually peaks at GeV energies in LBLs and at 
TeV energies in HBLs. 

The study of variability is one of the most powerful 
tools for revealing the nature of blazars and probing the 
various processes occurring in them. Based on their differ- 
ent timescales, the variability of blazars can be broadly di- 
vided into three classes, IDV, short-term variability (STV), 
and long-term variability (LTV). Variations in the flux of 
source up to a few tenths of mag nitude over a time scale 
of a day or less is known as IDV (|Wagner fc Witz"elll995l ) 
or microvariability, or intra-night optical variability. Varia- 
tions of days to a few months are often considered to be 
STV, while those from s everal months to many years are 
usually called LTV (e.g.. iGupta et al.l [20o3 l : both of these 
classes of variations for blazars typically exceed ~ 1 magni- 
tude and can exceed even 5 magnitudes. Over the last two 
decades, the optical variability of blazars has been exten- 



sively studied on diverse timescales (e.g., 


Heidt & Waenerl 


1996: Sillanpaa et al.lll996allbHGupta et al. 


2004 2008allblld, 


2009!: ICiprini et al.l 2003. 20071. and references therein). 



2 Rani et al. 



There are several theoretical models that might be able 
to explain the observed variability over wide time-scales 
for all bands, with the leading contenders all variants of 
mod els based upon shocks propagating down relativistic 
jets (iMarscher fc Gearll 19851; lOian et al. 199l|;lHughes et all 
ll99ll : lMarscher et all 19921 ; I Wagner fc Witzellll995l ). Some of 
the variability may arise from helical structures, precession 
or other geometrical effects occurr i ng within the jets (e.g.. 



Camenzind fc Krockenbergerlll992l ; iGopal-Krishna fc Wiital 
19921 ) and some of the radio variabilit y is due to extrinsic 



propagative effects (|Rickett et al.ll200 if ). Hot spots or other 
disturbances in or above accretion di sks surrounding the 
black holes at the centres of AGNs (e.g., lMangalam fc Wiital 
1 19931 ; IChakrabarti fc Wiital 119931 ) are likely to play a key 
role in the variability of non-blazar AGNs and might pro- 
vide seed fluctuations that could be advected into a rotating 
blazar jet and then be Doppler amplified. 

Despite the large amount of information we have about 
blazars that is very briefly summarized above, we still lack 
sufficient understanding of basic parameters of the emission 
regions, such as jet composition, a quantitative assessment 
of beaming parameters, or the processes leading to the origin 
of shocks in the jets. These physical quantities are obviously 
important in understanding the physics of jets and their 
emission regions and additional IDV studies leading to sta- 
tistically valid pictures of many blazars can help constrain 
them. 

In this paper we present the results of an extensive IDV 
studies of a sample of ten LBLs including six BL Lacs and 
four FSRQs. The work presented here is focused on intraday 
variations in the R passband magnitudes of these sources, 
which were the brightest blazars visible from ARIES, Naini- 
tal, India and PRL, Mount Abu India. We also compare 
our observational results to some of those presented in the 
literature. 

The paper is structured as follows. In Section 2, we 
present the observations and data reduction procedure. Sec- 
tion 3 provides our analysis and results. We present our dis- 
cussion and conclusions in Section 4. 



2 OBSERVATIONS AND DATA REDUCTION 

We have carried out optical R band photometric observa- 
tions of ten LBLs from September 2008 to June 2009 using 
two telescopes in India. The details of the telescopes and 
instruments used are given in Table 1 and the observation 
details are in Table 2. 

For doing image processing, or data pre-processing, we 
generated a master bias frame for each observing night by 
taking the median of all bias frames; the master bias frame 
was subtracted from all flat and source image frames taken 
on that night. Then the master flat in each passband was 
generated by median combine of fiat frames in that pass- 
band. Finally, the normalized master flat in each passband 
was generated. As usual, each source image frame was di- 
vided by the normalized master flat in the respective pass- 
band to remove pixel-to- pixel inhomogeneities (flat fielding). 
Finally cosmic ray removal was done from all source image 
frames. This pre-processing of the data was done by using 



standard routines in the Image Reduction and Analysis Fa- 
cility (IRAF) software. 

Our data analysis, or processing of the data, utilizes Do- 
minion Astronomical Observatory Photometr y (DAOPHOT 
II) so ftware to perform aperture photometry (|Stetsonlll987l . 
1992). We carried out aperture photometry with four differ- 
ent aperture radii, lxFWHM, 2xFWHM, 3xFWHM and 
4xFWHM. On comparing the results, we observed that 
aperture radii = 3xFWHM provided the best S/N ratio, 
and we have adopted that in this work. For all of the ten 
blazars, we observed more than three local standard stars. 
The magnitudes of the standard stars we used in the fields 
of our sources are given in Table 3. The multiple comparison 
stars were used to check that the usual standard stars were 
non-variable. We have used two non-varying standard stars 
from each blazar field and plotted their differential instru- 
mental magnitudes in Figs. 1—9. Finally, for the calibration 
of blazar data, we have used the one standard star that has 
a colour close to the blazar from those two standards stars. 
The calibrated light curves of the blazars are plotted in the 
same panel with the differential instrumental magnitudes of 
two standard stars. 



3 ANALYSIS AND RESULTS 
3.1 Variability Parameters 

We checked for the presence o f microvariability both by us- 
ing the F-test l)de D iego 2010) and the variability detection 
parameter, C, for the sake of comparison with earlier pa- 
pers, which nearly always used that approach. For two sam- 
ple variances, say Sq, for the quasar differential light curves 

b stS> 



F = 



(1) 



stS 



The F-statistic is compared with a critical value correspond- 
ing to the significance level set for the test. We have used 
the inbuilt F test code available in P0. A p-value of ^ 0.01 
(^99% significance level) is adopted for our variability de- 
tection criterion. 

The vari abilit y det ection parameter is defined as 
iRomero et alJ (e.g.. ll999T)) the average of Cland C2 where 



CI 



a(BL - star A) 
a(starA — starB) 



and C2 



a (BL - starB) 
(j (star A — starB) 



Here (BL — starA) and (BL — starB) are the differen- 
tial instrumental magnitudes of the blazar and standard 
star A and the blazar and standard star B, respectively, 
while <r(BL— starA) and <r(BL— starB) are the observational 
scatters of the differential instrumental magnitude of the 
blazar and star A and the blazar and star B, respectively. 
If C ^ 2.57, a conservative confidence level of a variability 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 R: A language and environment for statistical computing. R 
Foundation for Statistical Computing, Vienna, Austria. ISBN 3- 
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Optical Intraday Variability Studies of Ten Low Energy Peaked Blazars 3 



detection is > 99%, and we consider this to be a positive 
detection of a variation using this criterion. As noted by 
Ide Dieed (|20ld ) this C-statistic does not behave as a proper 
statistic should, but as it has been used in most of the IDV 
studies in literature, we also employed this test. 

The percentage variation in the intraday light curves of 
LBLs is calculated by usi ng the variability ampl itude pa- 
rameter, A, introduced bv lHeidt fc Wagner] (| 19961 ). defined 



A = 100 x y/(A max 



A r , 



2ff a (%), 



(3) 



where A max and A m in are the maximum and minimum mag- 
nitudes in the calibrated light curves of the blazar and a is 
the average measurement error of the blazar light curve. The 
calculated F-statistics, C- "statistics" and variability ampli- 
tude parameter, A, values are listed in Table 4. 

3.2 Structure Function 

The first order structure function (SF) is a very useful tool to 
search for periodi cities and timescales o f variability in time 
series data trains |Simonetti et al.l 1985). Here we give only 
a very brie f introduction to the method; for details refer to 
iRani et all (|2009h . The first order SF for a data train, a, is 
defined as 



Dl(k) 



S ^ w(i)w(i + k)[a(i + k) — a(i) 



(4) 



where k is the time lag, N a (k) = ^w(i)w(i + k), and the 
weighting factor, w(i) is 1 if a measurement exists for the 
i th interval, and otherwise. 

The behaviour of the first order SF can be simply sum- 
marized. The SF curves for AGN usually at first rise with 
time lag. Following this rising portion, the SF will then fall 
into one of the following classes: (i) if no plateau exists, any 
time scale of variability exceeds the length of the data train; 
(ii) if there are one or more plateaus, each one indicates a 
possible time scale of variability; and (iii) if that plateau is 
followed by a dip in the SF, the lag corresponding to the min- 
imum of that dip indicates a possible periodic cycle (unless 
such a dip is seen at a lag close to the maximum length of the 
data train, when it is probably an artifact). However, (iv) 
uncorrelated data produce a white noise behav iour, charac- 
terized by a constant slope (ICiprini et al.ll2003l ). 

We have carried out the SF analysis of all of those 
LCs which satisfy the variability detection criteria. Recently, 
some weaknesses of the SF method, including spurious indi- 
cat ions of timescales and periodici ties have been discussed 
by lEmmanoulopoulos et aL (2010), so we cross check the 
SF results by using the discrete correlation function (DCF) 
method. The timescales of variability calculated using the 
SF analysis are listed in Table 4. 



3.3 Discrete Correlation Functions 

The Discrete Correla tion function (DC F ) me thod was 
first introduced by lEdelson fc Krolikl (Il988l ) and it 
was later generalized to p rovide better error estimates 
ijHufnaeel fc Bregman|[l992l ). Here we give o nly a brief in- 
trodu ct ion to the method ; for details refer to lHovatta et al.l 
(5553), | Rani et all (|2009l ). and references therein. 



The first step is to calculate the unbin ned correlation 
(UDC F) using the given time series through (jHovatta et al.l 
120071 ) 



UDCF tJ = 



(a(i)-q)(6(i)-6) 



(5) 



Here a(i) and b(j) are the individual points in two time series 
a and b, respectively, a and 6 are respectively the means 
of the time series, and rj\ and <j\ are their variances. The 
correlation function is binned after calculation of the UDCF. 
The DCF method does not automatically define a bin size, 
so several values need to be tried. If the bin size is too big, 
useful information is lost, but if the bin size is too small, a 
spurious correlation can be found. Taking r as the centre of 
a time bin and n as the number of points in each bin, the 
DCF is found from the UDCF via 

DCF{r) = ^Y1 UDCF n( T )- 
The error for each bin can be calculated using 



(6) 



UDCF tJ - DCF(t) 



A DCF analysis is frequently used for finding the corre- 
lation and possible lags between multi-frequency AGN data 
where different data trains are used in the calculation (e.g. , 
IVillata et ai1l2004l : [Raiteri et all [20031 : iHovatta et al.l [20071 . 
and references therein). When the same data train is used, so 
a=b, there is obviously a peak at zero DCF, indicating that 
there is no time lag between the two, but any other strong 
peaks in the DCF give indications of variability timescales. 
The calculated t v values using DCF analyses are listed in 
Table 4. 



3.4 Intraday Variability (IDV) of Individual 
Blazars 

The R filter light curves (LCs) of the blazars in which 
significant variability has been detected, along with their 
corresponding SF and DCF analysis curves, are displayed 
in Figures 1—7; the remaining non-variable LCs of the 
blazars are plotted along with the corresponding curves 
of the standard stars used for comparison in Figures 8 
and 9. The complete observing log for the blazars in given 
in Table 3. The values of A, C, F-test along with any 
estimated timescales of variability found using SF and 
DCF analysis methods on the individual blazar LCs are 
listed in Table 4. A detailed multiband optical short term 
variability (STV) study of the fluxes and colours of all 
of these blaza rs over the s ame ti me period of observation 
is reported in IRani et"ail |2010b). There we showed that 
the colour versus brightness correlations seen in these 
sources support the hypothesis that BL Lacs tends to be 
bluer with increase in brightness while FSRQs shows the 
opposite trend. We now report some key individual results 
for each of our sources, placed in the context of earlier work. 

3C 66A: This is a low energy peaked blazar (LBL) at 
redshift z = 0.444 (lLanzetta et al.l Il993l ) and belongs to 
the class of BL Lac objects. Since its optical identification 
j| Wills fc Wills! 1 19741 ). the source has been regularly moni- 
tored at many observable frequencies, although less regularly 
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at radio fr e quencies (lAller et all Il992l ; iTakalo et all Il99fj) . 
iFan fc LirJ (|l999l . |2000) have studied the long-term optical 
and IR variability of the source and reported a variation of 
^1.5 mag at tim e scales of ~1 week to several years at those 
two frequencies. iBottcher et alj l|2005l ) reported a large mi- 
crovariability of ~0.2 mag within 6 hours; they also reported 
several major outbursts in the source separated by ~50 days 
and argued that the outbursts seem to have a quasi-periodic 
behaviour. At the end of 2007 the source was found to be in 
an optically active phase, which triggered a new Optical-IR- 
Radio Who le Earth Blazar Teles cope (WEBT) campaign on 
the source (jBottcher et alJl2009h . 

Our optical IDV observations of the source 3C 66A 
comprise a total of seven LCs, spanning a time period be- 
tween October 2008 and January 2009. During this period 
a ch ange of ~1 magni tude in brightness of the source is 
seen (R ani et al. I l2010bh . The source showed significant mi- 
crovariability only on October 22 and 26, 2008 (Fig. 1). 
There is a continuous fading trend of ~0.08 and ~0.06 mag- 
nitude, respectively on those two nights. The SF and DCF 
analysis of the LCs revealed that any timescale of variability 
in this source at those epochs is greater than the lengths of 
our observations. 

AO 0235+164: T he blazar AO 0235+164 at z = 0.94 
llNilsson et all 1 19961) was classified as a BL Lac object by 
ISpinrad fc Smithl l| 19751 ). Over the past few decades this 



and at all frequencies 
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raiaDerumallll995!: 


iHeidt & Wagnerlll99e 


iRaiteri et al.1 2001 


) and a very high 



fractional polarization of ~40% has been rep orted in the 
sourc e both at visible and IR frequencies (jlmpev et al.l 
Il982h . B y analysing 25 years (1975-2000) of optical and ra- 
dio data. IRaiteri et aD ( 200 ll ) argued that the source seemed 
to have an optical outburst period of ~5.70 years but the 
expected outburst in 2004 was not detected b y a 2003-2005 
multiw avelength WEBT observing campaign ()Raiteri et al.l 
2006a). A more detailed long term optical data analysis sug- 
gested a possible out burst period of ~8 years in this source 
l)Raiteri et al.ll2006bl ) and this period was supported by sub- 
sequent observations iGupta et ail ( |2008<J > . Strong IDV flux 
variations of 9.5% and 13.7% during tw o nights were ob - 
served bv lGupta et~aD <|2008<j + Recently. iRani et all ()2009l ) 
reported the possible presence of nearly periodic fluctua- 
tions, with a timescale of ~17 days, in a 12.7 year long X-ray 
light curve of AO 0235+164 obtained by the All Sky Moni- 
tor (ASM) instrument on the Rossi X-ray Timing Explorer 
(RXTE) satellite. 

We observed three IDV LCs of the source AO 0235+164 
between October 2008 and January 2009. The brightness of 
the source decayed by ~2.2 magnitude during this period 
l)Rani et al.ll2010bh . We found a significant flux variation in 
two out of three LCs. A continuous fading trend of ~0.13 
magnitude and both brightening and decaying of ~0.04 mag- 
nitude were observed on October 20 (Fig. 1) and October 23, 
2008 (Fig. 2) , respectively. A possible timescale of variability 
is ~5.2 hr (from the SF analysis) for the LC observed on 20 
October, while any such timescale exceeds the length of our 
observation on 23 October. However this putative timescale 
is not supported by the DCF analysis. 

PKS 0420-014: The blazar PKS 0420-014 is classified as 
a FSRQ and has a redshift of 0.915. It has been observed 
in optical bands since 1969. Several papers have reported 



multiple optically active and bright phases of the source 
and perhaps regular ma j or flaring cycles (e. g., IVillata et al.l 
1 19971 : 1 Webb et al] 1 19981: IRaiteri et all 1 19981 . and references 
therein). IWebb et al.l (119981 ) reported that there were in- 
creases of ~ 2—3 magnitudes during the active phases of 
this blazar during their observa tions that stretched fro m De- 
cember 1969 to January 1986. IClements et"afl (|l995l ) have 
reported variations of Amag = 2.8 mag with a time scale of 
~22 years. 

We found a variation of >0.1 magnitude within 2 hrs 
in the brief optical LCs of the source on both of the days 
of observation in October and December 2008. During this 
perio d the source brigh tened by a factor of ~0.7 magni- 
tude (|Rani et alj|2010br ). The nominal timescale of variabil- 
ity (from the peak of the SF) is 0.12 hrs for the LC observed 
on December 26 and multiple dips might indicate a possible 
periodicity around 0.18 hrs, which is weakly supported by 
the DCF curve. The other LC is irregular and shows no hint 
of a timescale (Fig. 2). 

S5 0716+714: The blazar S5 07 1 6+714 is classified 
as a BL Lac object. iNilsson et all (120081 ) made a re- 
cent claim of redshift determination of the source to be 
z — 0.31 ± 0.08. This source has been extensively stud- 
ied at all observable wavelengths from radio to 7-rays on 
diverse time scales (|Wagner et al. 1 1 19901; IHeidt fc Wagner) 



1996; Villata et al. 2000; Raitcri et al. 2003; Montagni et al 



2006; Fosc hini et all 2006; Os torero et a l. 2006; Gupta et al 



2008al ldV This source is one of the brightest BL Lacs in op- 



tical bands with an IDV duty cycle of nearly 1. Unsurpris- 
ingly, it has been the subject of s everal optical monitor- 
ing campaigns on IDV timescales jHeidt fc Wagner 1996; 
iMontagni et al.ll2006l ; [Gupta et al.ll2008d). This source has 
shown five major optical outbursts ( Gupta et al.ll2008d ) at 
intervals of ~3.0±0.3 years. High optical polar izations of ~ 
20% - 29% has also been observed in the so urce |Takalo et al.l 
ll994l ; lFan et al.1 19971 ). iGupta et al.l j|2009l ) reported good ev- 
idence for nearly periodic oscillations in a fe w of the intraday 
optica l light curves of the source observed bv lMontagni et all 
(2006). Good evidence of presence of a ^15 minute peri- 
odic oscilla t ion at optical frequencies has been reported by 
IRani et al.l (|2010al ). 

Our optical IDV observations of S5 0716+714 spans a 
time period from October 2008 to January 2009. The source 
brightened by a fa ctor of ^2 magnitude during this period 
l|Rani et al.ll2010bl ). We found significant microvariability of 
~0.1 magnitude in three out of four LCs of the source. The 
LCs observed on December 24, 2008 and January 03, 2009 
show continuous fading trends trend of the order of ~0.1 
magnitude, though the former is abrupt while the latter 
is gradual. Both fading and brightening and fading trends 
of ~0.05 magnitude were observed over just a few minutes 
on December 23, 2008. The calculated possible variability 
timescales are listed in Table 4, but the lack of agreement 
between the SF and DCF possibilities leads us to discount 
their reality. 

PKS 0735+178: The blazar PKS 0735+718 has b een clas- 
sified as a BL Lac object l|Carsw cll et al. 197 4]). Papers 
concerning its redshift determination ( Carswell et al.l Tl9T4l : 
iFalomo fc Ulrichll2000l) had set a lower limit of z > 0.4 and 
z — 0.424 was reported for the so urce using a HST snap- 
shot image (|Sbarufatti et al]|2005h . Since its optical iden- 
tification, the source has been extensively observed across 
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the whole elect r omagnetic spectrum llTerasranta et ai]|2004l ; 
iGu et al.1 [200^ : iGupta et al.l |2008(J : ICiprini et al.l 1200^ . A 

periodicity of ~14 years has been suggested to be present 
in the sour ce using a centu ry long, but still sparse, optical 
light curve l|Fan et al.ll 19971 ). Optical variability o n IDV and 
STV t imescales has been observed for 0735+ 1 78 dXie et al 



1992) ; iMassaro et al.1 Il995l; iFan et al.l Il997l ; IZhang et a". 
2004 ; ICiprini et al.ll2007l ; IGupta et al.l l2008d l . A significant 



amount of polarization (~ 1% to 30%) ha s been observed in 
the source bo t h at optical and IR bands dMead et al.l 199C ; 



Takalol Il99ll; iTakalo et alj Il992bl : IValtaoia et~ 
19931 ; iTommasi et al.ll200ll ) 



1991a 



Our IDV observations of the source PKS 0735+718 
comprise four LCs taken between December 2008 and Jan- 
uary 2009. The brightness of the source chang es by ~0.6 
magnitude during this period (jRani et al.ll2010q ). We found 
significant microvariations of an order of ~0.1 magnitude in 
three out of four observed LCs of the source (Fig. 4). The 
only conceivable hints of timescales for variability from the 
SFs are ~2.3 hrs for the LC observed on December 28, 2008 
and ~0.58 hrs for January 04, 2009. However, since both of 
these peaks are close to the total lengths of the observations 
they are not likely to be real, nor supported by DCF results. 
OJ 287: The blazar OJ 287 at z = 0.306 is one of 
the most extensively observed and best studied BL Lac 
objects. It is also among the very few AGN's for which 
more than a century o f optical observations are available 
llSillanpaa et al l Il996allbl; IFan et al.l 1 19981 ; lAbrahaml |2000| ; 
IGupta et al.ll2008d; IFan et alJ^OOgf T Using the binary black 
hole model IjSillanpaa, et al.l 1 1988 ) for the long-term opti- 
cal light curve of the source, an outburst with a predicted 
~12 year p eriod was detected in the source by the OJ -94 
programme (|Sillanpaa et al. l ll996al ; IValtonen et al.ll200sh . A 
very high optical polarization and variability in the degree 
and angle of polarizat ion has been also reported for OJ 
287 (Hi mov et al]|2002r i . The observational properties of the 
so urce from radio to X -ray energ y bands have bee n reviewed 
bv lTakalo et~aH (|l994T ). Recently. IFan et al.1 (|2009h , reported 
large variations in the source of AV = 1.96 mag, AR = 2.36 
mag, and Al = 1.95 mag during their observations spanning 
2002 to 2007. 

Our observations of OJ 287 span a period from De- 
cember 2008 to January 2009, and during this time the 
brightn ess of the source c hanged significantly by ~1 mag- 
nitude dRani et al.1 l2010bl ). The source showed significant 
microvariations only in one out of three observed nightly 
LCs during which the brightness of the source continuously 
faded by ~0.08 magnitude within 2 hrs (Fig. 5). The SF 
and DCF analysis showed that any timescale of variability 
is longer than the timescale of observation. 
3C 273: The FSRQ 3C 273 was th e first quasar dis- 
covered and has a redshif t of 0.158 dSchmidtl 1 1963ft . It 
is categorized as a LBL dNieppola et al.l I2006T ) and its 
spectral energy distribution, correlations among differ- 
ent energy band flares and the approaching jet orien- 
tation have been extensively studied at all wavelengths. 
There are many papers covering the obse rvational prop- 
erties of the source in the optical band (lAngiond 
Sitko et al.lll982l;ICorso et al.|[l985l. Il986l: iMqles et al.l 



1971 



1986 



Hamuv fe Mazal 19871; ISillanpaa, et al.lll99ll; IValtaoia et al 
1991b; Takalo ct al. 1992a b; Elvis ct al. 1994; Lichti ct al 



19951 ; iGhosh et al.1 120001 ') . An analysis of the optical light 



curve of 3C 273 spanning over 100 years c an be interprete d 
to suggest a LTV timescale of ~13.5 years (|Fan et al.ll200ll ). 
Recently, the short-term optical variability and colour in- 
dex pr operties of the source have been studied bvlDai et al.l 
l|2009ft . 

Our IDV observations of the source 3C 273 span the pe- 
riod from December 2008 to April 2009 during which a total 
of eight LCs were obtained. There is no change in overa ll 
flux of the source during this period l|Rani et al.l l2010bft . 
This source showed microvariations (of ~0.05 magnitude) 
in only two out of eight observed LCs (Fig. 5). On the night 
of 19 April 2009 there is a hint of a timescale of variability 
of ~6 hrs from both the SF and DCF approaches. 
PKS 1510-089: The blazar PKS 1510-089 is classified as 
a FSRQ and has z — 0.361. It also belongs to the category of 
highly polarized quasars. Significant optical fl ux variations 
in the source were first reported by iLul l|l972ft over a time 
span of ~5 years. The historical light curve of the source 
shows a large variation of AB = 5.4 mag during an out- 
burst in 1948 after which it faded by ~2.2 mag within 9 
days (|Liller fc Lillerlll975l ). Strong variations on IDV time 
scales also have been report ed for PKS 151 0—089; e.g., AR 
= 0.65 mag within 13 min. |Xie et al.ll200lft , AR = 2.0 mag 
in 42 min jDai et al.l l2001h . AV = 1.68 mag in 60 min. 
(|Xie et al.ll2002aft . In the optical light curves of this source, 
a few deep minima have been observed on dif ferent days 
jXie et al.l200ll ; lDai et al.l200ll ; lxie et al.ll2002bft . that nom- 
inally correspond to a time scale of ~42 min, though no more 
than 3 such dips were ever seen in a single night. Nonethe- 
less, an eclipsing binary black hole model was act ually pro- 
posed to explain the occurrences of these minima (|Wu et al.l 
12005ft . Other observations by this group have yielded a claim 
of another possi ble time scale between minima of ~89 min 
|Xie et al.ll2004T ). 

We carried out IDV observations of the source from 
April to June 2009. There is a large change in the bright- 
ness of the source during this period , with AR ma9 = 1.5 
ijRani et al.ll2010bft . We found significant microvariations of 
~0. 05-0. 08 magnitude in three out of four LCs (Fig. 6). We 
found that any timescale of variability is larger than the 
timescale of observations, except perhaps for the LC ob- 
served on April 19, 2009 for which it is formally ~0.6 hrs 
from the SF curve and ~00.5 hi from the DCF; however, this 
LC has too few points to allow the production of a crisp SF 

or DCF, so this evidence is very weak. 

BL Lac: The object BL Lac at z = 0.069 (|Miller et al.lll978l ) 
is the archetype of its class. Observations over the past few 
decades have showed that its optical and radio emissions are 
highly variable and polarized and the polarization at those 
widely separated freque ncies is found to be strongly corre- 
lated l|Sitko et al l 1 1985ft . It is among the very few sources 
for which more t han 100 yea r s of optical data i s available 
in th e literature (|Shenl Il970l ; I Webb > et all QUI; IFan et al.1 
1998). An optical variation of AB = 5.3 mag and a possible 
periodicity of ^ 14 years has b een reported for BL L ac by 
IFan et al.1 1| 19981 ). Very recently. iNieppola et all |2009l ) have 
studied the long term variability of the source at radio fre- 
quencies and generalized the shock model that can explain 
it. 

Our IDV observations of the source BL Lac were made 
between September 2008 and June 2009. The brightness 
of the source faded by ~1.6 magnitude during this pe- 
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riod (|Rani et al.ll2010bl ). BL Lac faded by ~0.05 magnitude 
within 1 hr of observation on September 04, 2008 which had 
a nominal variability timescale of ~2.5 hrs according to both 
the SF and DCF plots (Fig. 7). A continuous rising trend 
of ~0.1 magnitude in the LC of the source was observed on 
June 21, 2009, and the calculated timescale of variability for 
the LC that night exceeds the time period of the observa- 
tions. 

3C 454.3: A FSRQ at a redshift of 0.859, 3C 454.3 
is among the most intense and variable sources. The 
source has been detected in the flaring state in July 2007 
and July 2008 at 7-ray frequencies and those flares have 
been found to be well correlated with optical and longer 
wavelength flares (fGnisellini et al.ll2007l ; iRaiteri et ai1l2008l ; 
IVillata et alj 12007 ). The long term observational proper- 
ties of the source at optical and radio frequencies have 
been well studied through multiwavelength campaigns (e.g., 
IVillata et al.ll2006l, \20m. The IDV of the source was re- 
cently studied by Gupta et al.l (|2008d ). They have reported 
that the amplitude varied by ~5 - 17% during their ob- 
serving span. An extraordinary flaring activity above 100 
MeV has been rep orted in the source in December 2009 
|Striani et al.ll2010l ). 

We observed two IDV LCs of the source on October 24 
and 28, 2008. During the peri od the brightness of the source 
increases by ~0.4 magnitude (|Rani et al ]|2010bl ). The source 
showed significant microvariations of ~0.13 magnitude on 
October 24, 2008 (Fig. 7) while no significant variations has 
been detected for the LC observed on October 28, 2008. 
The SF and DCF analysis revealed that any timescale of 
variability exceeds the length of the observation on the first 
of those nights. 



4 DISCUSSION AND CONCLUSIONS 

We have carried out optical R band IDV observations of ten 
LBLs spanning over a time period of 2008 September to 2009 
June. The sources PKS 0420-014 and PKS 1510-089 were 
in faint states; AO 0235+164, BL Lac and 3C 454.3 were 
possibly in post-outburst states; S5 0716+714 and 3C 66A 
were in pre-outburst states, while PKS 0735+178 and OJ 
287 were in some inter mediate states during the observing 
run l|Rani et al.l l2010bl ) . In our search of microvariations in 
ten LBLs we found significant IDV in 21 out of 40 observed 
LCs; so the calculated duty cycle of IDV in LBLs during 
our observing run is ~52%. We performed the SF and DCF 
analysis to calculate the nominal timescales of variability; 
however we found that in most of the cases this timescale of 
variability is longer than the length of observations and in 
a substantial majority of the cases where the SF indicated 
a possible timescale the DCF did not support it. 

The blazar emission mechanism in the outburst state is 
quantitatively understood by relativistic shocks propagating 
through a relativistic jet of plasma. In general, blazar emis- 
sion in the outburst state is non-thermal Doppler-boosted 
emission fr om jets enhanced by that arising from shocks in 
the flows. llBlandford fc Reeslll978l ; IMarscher fe Gear! 1 19851 ; 
iMarscher et al. 19921 ). The other models of AGNs that can 
explain the IDV in any type of AGN are optical flares, dis- 
turbances or hot spots on the a ccretion disk surrounding 
the black hole of the AGN (e.g., iMangalam fc Wiital fl993l 



and references therein). Models based on the instabilities 
on the accretion disk could convinancibly yield blazar IDV 
only when the blazar is in the very low state. When a blazar 
is in the low state, any contribution from the jets, if at all 
present, is very weak. So, we consider that the observed 
IDV in the sources AO 0235+164, BL Lac, 3C 454.3, S5 
0716+714, 3C 66A, 0735+178 and OJ 287 is almost cer- 
tai nly related to a shock propagating through a relat ivistic 
jet (|Blandford fc Konigl|[l^79l ; IMarscher fc Gearll 19851 ). Tur- 
bulence behind a shock propagating down s uch a jet is a ver y 
feasible way to explain the observed IDV (|Marscherlll996l ). 
Since the blazars PKS 0420-014 and PKS 1510-089 were 
observed in relatively faint states, there is a chance that the 
observed optical IDV in these source may be because of hot 
spots or any other enhanc ed emission on the accretion disk 
|Mangalam fc Wiital fl993l ). 

In one source, 3C 273, we are unable to classify the state 
of the source since this blazar has been in an essentially 
steady state for last six years |Dai et alj 120091 ; iRani et"all 
l2010bl ). This source showed significant microvariations in 
two out of eight observed LCs and on both of days of obser- 
vation the brightness of the source follows both rising and 
fading trends of ~0.05 magnitude. Whatever may be the 
mechanism responsible for the origin of of microvariability 
in this source, it does not seem to be strong enough to in- 
troduce day-to-day variations in the flux of the source. 

It is worth noting that an extrinsic mechanism can also 
be responsible for some of the observed IDV in blazars. Ex- 
trinsic IDV could be caused by refractive interstellar scin- 
tillation, which is only relevant in low-frequency radio ob- 
servations, or microlensing, which is achromatic. We note 
that the blazar AO 0235+ 164 has two foreground galaxies 
at z = 0.524 and z = 0.8 51 (|Cohen et al.ll 19871 ; iNilsson et all 
1 19961 ; IWebb et al] I2OO0I ). The flux of this source is con- 
taminated and absorbed by foreground galaxies, the stars 
of which can act as gravitational microlenses. Thus the 
observed optical IDV in AO 0235+164 could arise, at 
least partially, from gravitational microlensing. Using two 
separated telescopes to s i multa neously observe 0716+714 
iPollock. Webb fc Azarnial l|2007h showed that instrumental 
and atmospheric effects cannot account for the microvaria- 
tions they measured for that blazar. 

We found significant IDV in 21 out of 40 observed LCs 
of ten LBLs; so the calculated duty cycle of IDV in LBLs 
during our observing run is ~52%. The average duration of 
our observation was 3.7 hrs per LC. The SF and DCF curves 
revealed that in ~60% cases any timescale of variability is 
longer than the timescale of observations. In quite a few 
cases there were hints of timescales in the data from the SF 
plots, but only a few of those hints were supported by the 
DCF plots. 

Extensive IDV studies of different subclasses of AGNs 
revealed that the occurrence of IDV in blazars observed on 
a timescale of <6 hrs is ~60— 65% and if the blazar is ob- 
served more tha n 6 hrs then the possibility of ID V detec- 
tion is 80-85% (|Carinilll990l ; iGupta fc Joshill2005l . and ref- 
erences therein). If we consider observations over days to 
months, i.e., at STV timescales then the observed duty cy- 
cle of variations is >92% (e.g. IRani et al.1 l2010bl ) and at 
LTV timescales it is almost 100%, which confirms that the 
probability of detection of variability in blazars rises with 
the duration of the observations. Although the duty cycle 
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from our current obse rvations is less than that reported by 
iGuptafc Jo"shl(|2005l ), since the average duration of our ob- 
servations is <4 hrs, this is unsurprising. 
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Table 1. Properties of Telescopes and Instruments 



Site: 


ARIES Nainital 


PRL Mount Abu 


Telescope: 


1.04-m RC Casscgrain 


1.20 m Casscgrain 


CCD model: 


Wright 2K CCD 


Andor EMCCD 


Chip size: 


2048 X 2048 pixels 


2048 x 2048 pixels 


Pixel size: 


24 X 24 fim 


13 x 13 ;um 


Scale: 


0.37" /pixel 


0.17" /pixel 


Field: 


13' x 13' 


3' x 3' 


Gain: 


10 e-/ADU 


5 e-/ADU 


Read Out Noise: 


5.3 e~ rms 


4.9 e~ rms 


Binning used: 


2x2 


2x2 


Typical seeing : 


l"to 2.8" 


l"to 2.6" 



Table 2. Observation Log 



Blazar Name 


Date of Observation 


Telescope 


Filter 


Data Points 


Duration (hrs) 


3C 66A 


2008 Oct 22 


A 


R 


68 


3.86 




2008 Oct 26 


A 


R 


72 


3.50 




2008 Dec 23 


B 


R 


90 


1.45 




2008 Dec 24 


B 


R 


113 


1.85 




2008 Dec 27 


B 


R 


417 


3.46 




zUUo Dec ZO 


B 


R 


203 


2.23 




2009 Jan 03 


B 


R 


320 


3.54 


AO 0235+164 


2008 Oct 20 


A 


K 


98 


6.30 




2008 Oct 23 


A 


R 


45 


2.60 




2008 Dec 26 


B 


R 


248 


1.23 


PKS 0420—014 


2008 Oct 23 


A 


R 


18 


1.94 




2008 Dec 26 


B 


R 


29 


0.73 


S5 0716+714 


2008 Oct 24 


A 


K 


25 


1.35 




2008 Dec 23 


B 


R 


114 


0.40 




2008 Dec 24 


B 


R 


292 


1.62 




2009 Jan 03 


B 


R 


2685 


3.73 


PKS 0735+178 


2008 Dec 23 


B 


R 


90 


0.63 




2008 Dec 28 


B 


R 


300 


3.20 




2009 Jan 04 


B 


R 


50 


1.02 




2009 Jan 20 


A 


R 


67 


3.98 


OJ 287 


2008 Dec 26 


B 


R 


70 


0.57 




2009 Jan 20 


A 


R 


60 


3.87 




2009 Jan 22 


A 


R 


70 


3.35 


3C 273 


2008 Dec 23 


B 


R 


130 


0.72 




2009 Jan 22 


A 


R 


90 


3.88 




2009 Feb 25 


A 


K 


101 


3.61 




2009 Mar 24 


A 


K 


36 


1.38 




2009 Apr 01 


A 


K 


110 


5.60 




2009 Apr 18 


A 


R 


91 


3.54 




2009 Apr 19 


A 


K 


195 


7.42 




2009 Apr 27 


A 


R 


66 


3.31 


PKS 1510-089 


2009 Apr 17 


A 


R 


77 


3.82 




2009 Apr 19 


A 


R 


22 


1.09 




2009 Apr 27 


A 


R 


64 


4.47 




2009 June 21 


A 


R 


68 


4.41 


BL Lac 


2008 Sep 04 


A 


K 


83 


5.10 




2008 Oct 26 


A 


R 


73 


4.25 




2009 June 21 


A 


R 


62 


2.98 


3C 454.3 


2008 Oct 24 


A 


R 


55 


3.01 




2008 Oct 28 


A 


R 


65 


4.25 



A : 1.04 m Sampuranand Telescope, ARIES, Nainital, India 
B : 1.20 m Telescope, PRL, Mount Abu, India 
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Table 3. Standard Stars in the Blazar Fields 



Source 


Standard 


R magnitude 


Refrences 


Name 


star 


(error) 




3C 66A 


1 


13.36(0.01) 


5 




2 


14.28(0.04) 


5 




3 


15.46(0.12) 


5 




4 


12.70(0.04) 


6 




5 


13.62(0.05) 


6 


AO 0235+164 


1 


12.69(0.02) 


1 




2 


12.23(0.02) 


1 




3 


12.48(0.03) 


1 




6 


13.64(0.04) 


6 




8 


15.79(0.10) 


1 




CI 


14.23(0.05) 


6 


PKS 0420—014 


1 


12.09(0.03) 


1 




2 


12.80(0.02) 


5 




3 


12.89(0.01) 


5 




4 


14.47(0.01) 


5 




5 


14.37(0.03) 


1 




6 


14.70(0.03) 


1 




7 


14.91(0.03) 


1 




8 


15.46(0.03) 


4 




9 


15.58(0.04) 


1 


So 0716+714 


1 


10.63(0.01) 


2 




2 


11.12(0.01) 


2, 3 




3 


12.06(0.01) 


2, 3 




4 


12.89(0.01) 


2 




5 


13. 18(0. 01) 


2, 3 




6 


13.26(0.01) 


2, 3 




7 


13.32(0.01) 


2 




8 


13.79(0.02) 


2 


PKS 0735+178 


A 


13.14(0.05) 


1 




C 


13.87(0.06) 


1 




D 


15.45(0.06) 


1 


OJ 287 


2 


12.46(0.05) 


1 




4 


13.72(0.06) 


1 




10 


14.26(0.06) 


1 




11 


14.67(0.07) 


1 


3C 273 


C 


11.30(0.04) 


1 




D 


12.31(0.04) 


1 




E 


12.27(0.05) 


1 




G 


13.16(0.05) 


1 


PKS 1510—089 


1 


11.23(0.03) 


5 




2 


12.95(0.03) 


5 




3 


13.98(0.09) 


5 




4 


14.34(0.05) 


5 




5 


14.35(0.05) 


1 




6 


14.61(0.02) 


1 


BL Lac 


B 


11.93(0.05) 


1 




C 


13.69(0.03) 


1 




H 


13.60(0.03) 


1 




K 


14.88(0.05) 


1 


3C 454.3 


A 


15.32(0.09) 


6, 9 




B 


14.73(0.05) 


6, 9 




C 


13.98(0.02) 


9, 10 




D 


13.22(0.01) 


9, 10 




E 


14.92(0.08) 


6, 9 




F 


14.83(0.03) 


4, 9 




G 


14.83(0.02) 


4, 9 




H 


13.10(0.04) 


6, 9 




CI 


15.27(0.06) 


6 



a l. l lSmith et al.lll985l); 2 llVillata et alJll99Sh : 3. l lGhisellini et al.lll997l) ; 4. llRaiteri et al.lll99Sf) ; 5. llSmith fc BaIoneklll998T ); 6. 
ilFiorucci &: Tost3ll996h ; 7. Craine E.R.:Handbook of Quasistellar and BL Lacertae Objects. lAligione!Il97ll ); 8. 
http://www.lsw.uni-heidelberg.de/projects/extragalactic/charts/2251+158.html 
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Table 4. IDV Results 



Blazar Name 


Date of 


C - Test 


F - 


Test 


V 


A (%) 








Observation 


value 


F-value 


p- value 






SF (hrs) 


DCF (hrs) 


3C 66A 


2008 Oct 22 


3.59 


3.65 


3.1c" 7 


V 


8.1 


>3.86 


>3.86 




2008 Oct 26 


2.65 


7.03 


1.9e" 14 


V 


5.7 


>3.50 


>3.50 




2008 Dec 23 


0.78 


0.41 


0.06 


NV 










2008 Dec 24 


3.31 


0.61 


0.02 


NV 










2008 Dec 27 


0.88 


0.58 


0.26 


NV 










2008 Dec 28 


0.75 


0.37 


0.3 


NV 










2009 Jan 03 


0.83 


0.63 


0.35 


NV 








AO 0235+164 


2008 Oct 20 


8.32 


76.84 


2.2e~ 16 


V 


13.4 


5.20? 


>6.30 




2008 Oct 23 


3.23 


10.53 


1.4e- 14 


V 


4.2 


>2.60 


>2.60 




2008 Dec 26 


1.09 


1.43 


0.04 


NV 








PKS 0420-014 


2008 Oct 23 


5.44 


31.44 


6.6e -15 


V 


14.4 








2008 Dec 26 


5.42 


31.24 


6.1e- 15 


V 


12.3 


0.12,0.18 


0.18 


S5 0716+714 


2008 Oct 24 


0.92 


1.57 


0.27 


NV 










2008 Dec 23 


3.47 


7.77 


2.2c~ 16 


V 


9.1 


0.0014 


0.0028 




2008 Dec 24 


4.51 


3.07 


2.2c- 16 


V 


14.6 


>1.62 


>1.62 




2009 Jan 03 


3.71 


1.26 


7.7c- 7 


V 


31.6 


3.28 


>3.73 


PKS 0735+178 


2008 Dec 23 


3.14 


9.49 


2.2c- 16 


V 


11.1 








2008 Dec 28 


4.22 


18.36 


2.2c- 16 


V 


19.3 


2.90 


2.30 




2009 Jan 04 


2.31 


4.64 


3.1e" 7 


V 


9.8 


0.58 


0.60 




2009 Jan 20 


1.30 


1.99 


0.055 


NV 








OJ 287 


2008 Dec 26 


1.11 


1.48 


0.10 


NV 










2009 Jan 20 


0.88 


1.20 


0.47 


NV 










2009 Jan 22 


3.21 


11.13 


2.2e~ 16 


V 


7.9 


>3.35 


>3.35 


3C 273 


2008 Dec 23 


0.91 


0.73 


0.08 


NV 










2009 Jan 22 


0.81 


0.65 


0.05 


NV 










2009 Feb 25 


0.75 


0.49 


0.06 


NV 










2009 Mar 24 


0.76 


0.50 


0.05 


NV 










2009 Apr 01 


2.59 


2.58 


1.2c- 6 


V 


5.2 


>5.2 


>5.60 




2009 Apr 18 


1.59 


0.25 


0.11 


NV 










2009 Apr 19 


3.74 


2.22 


4.7e- 8 


V 


4.5 


6.0 


6.0 




2009 Apr 27 


0.86 


0.84 


0.49 


NV 








PKS 1510-089 


2009 Apr 17 


2.67 


2.31 


0.0003 


V 


9.3 


>3.82 


>3.82 




2009 Apr 19 


2.72 


7.11 


3.3e- 5 


V 


7.6 


0.60 


0.50 




2009 Apr 27 


1.27 


1.54 


0.09 


NV 










2009 June 21 


4.46 


2.65 


0.0001 


V 


8.3 


>4.41 


>4.41 


BL Lac 


2008 Sep 04 


2.73 


5.71 


l.Oe- 13 


V 


8.7 


2.5 


2.5 




2008 Oct 26 


1.05 


1.56 


0.06 


NV 










2009 June 21 


2.87 


7.34 


4.9e" 13 


V 


12.9 


>2.98 


>2.98 


3C 454.3 


2008 Oct 24 


5.14 


22.95 


2.2c- 16 


V 


16.3 


>3.01 


>3.01 




2008 Oct 28 


1.42 


2.09 


0.01 


NV 
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Figure 1. R band optical IDV LCs of the blazars 3C 66A and AO 0235+164 and their respective SFs and DCFs. 
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Figure 2. R band optical IDV LCs of the blazars AO 0235+164 and PKS 0420-014 and their respective SFs and DCFs. 
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Figure 3. R band optical IDV LCs of the blazar S5 0716+714 and their respective SFs and DCFs. 
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Figure 4. R band optical IDV LCs of the blazar PKS 0735+178 and their respective SFs and DCFs. 
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Figure 5. R band optical IDV LCs of the blazars 3C 273 and OJ 287 and their respective SFs and DCFs. 
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Figure 6. R band optical IDV LCs of blazar PKS 1510-089 and their respective SFs and DCFs. 
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Figure 7. R band optical IDV LC of the blazars BL Lac and 3C 454.3 and their respective SFs and DCFs. 
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Figure 8. The R filter LCs of several blazars during which no significant variability has been detected, plotted along with the differential 
magnitudes of standard stars. 
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Figure 9. The R filter LCs of several blazars during which no significant variability has been detected, plotted along with the differential 
magnitudes of standard stars. 



